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The  purpose  of  this  paper  is  to  show  that  any  generalized  network 
problem  whose  matrix  does  not  have  full  row  rank  is  equivalent  to  a 
pure  network  problem.  Furthermore,  wc  develop  an  efficient  pro¬ 
cedure  which  may  be  used  simultaneously  to  determine  the  rank  of 
the  incidence  matrix  anti  transform  (if  possible)  the  generalized  net¬ 
work  problem  into  a  pure  problem.  Transforming  such  problems  into 
pure  network  problems  provides  several  advantages,  such  as: 

1)  increasing  computational  efficiency 

2)  eliminating  roundoff  error 

2 )  providing  a  simple  check  for  infeasibility 

4)  extending  the  results  in  the  literature  on  pure  networks  to 
appropriate  generalized  networks. 
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1 . 0  Introduction  and  Conclusion 

The  purpose  of  this  paper  is  to  show  that  any  generalized  network 
problem  whose  incidence  matrix  does  not  have  full  row  rank  is  equivalent 
to  a  pure  network  problem.  Furthermore,  we  develop  an  efficient  pro¬ 
cedure  which  may  be  used  simultaneously  to  determine  the  rank  of  the 
incidence  matrix  and  transform  (if  possible)  the  generalized  network 
problem  into  a  pure  network  problem. 

An  important  generalized  network  problem  that  can  be  transformed 
into  a  pure  network  problem  by  our  procedure  is  the  stock  cutting  problem 
of  mathematical  programming.  In  general,  a  number  of  generalized  net¬ 
work  problems  have  a  structure  that  permits  them  to  undergo  this  trans¬ 
formation.  For  instance,  upon  organizing  the  incidence  matrix  of  a  generalized 
transportation  problem  in  echelon-diagonal  form,  any  generalized  transpor¬ 
tation  problem  with  repeating  echelons  and  repeating  diagonals  has  the  appro¬ 
priate  rank. 

Transforming  generalized  network  problems  into  pure  network  problems 
provides  obvious  computational  advantages.  While  the  generalized  network 
problem  is  next  in  order  after  pure  networks  or  distribution  problems  with 
respect  to  ease  of  computation,  the  .jump  in  degree  of  difficulty  is  such  that 
there  exist  no  algorithms  for  the  former  class  of  problems  comparable  in 
speed  or  efficiency  to  those  for  the  latter.  In  fact,  no  computation  times 
for  special  purpose  algorithms  for  solving  generalized  network  and 
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transportation  problems  (1,2,6,  11,  12,  13]  have  ever  been  published.  More¬ 
over,  none  of  the  major  computer  manufacturers  have  indicated  the  availa¬ 
bility  of  special  purpose*  codes  for  solving  generalized  network  and  trans¬ 
portation  problems  in  their  libraries  or  commercial  packages.  Thus  we 
conclude  that  the  existence  of  such  codes  must  be  rare*,  and  for  all  practical 
purposes  the  only  way  currently  available  for  solving  a  generalized  network 
or  transportation  problem  is  via  a  general  purpose  linear  programming 
code.  This  situation  magnifies  the  discrepancy  in  computational  effort 
required  to  solve  pure  and  generalized  network  problems.  As  reported  in 
[9],  transportation  problems  can  be  solved  at  least  30  times  faster  by  a 
special  purpose  code  (utilizing  the  techniques  of  [7,  8])  than  by  general 
purpose  linear  programming  methods  (utilizing  sophisticated  procedures  for 
exploiting  sparse  matrices).  This  difference  would  undoubtedly  be  accen¬ 
tuated  when  the  general  purpose  code  is  applied  to  a  generalized  network 
structure. 

There  are,  of  course,  other  advantages  to  transforming  a  generalized 
problem  into  a  pure  network  problem  when  this  is  possible.  First,  by  using 
appropriate  specialized  algorithms,  no  round-off  error  is  introduced  while 
solving  the  pure  network  problem.  Second,  if  the  total  supply  of  the  pure 
network  problem  is  not  equal  to  the  lotal  demand,  then  one  can  conclude  without 
further  effort  that  the  original  problem  lacks  a  feasible  solution.  (Similarly, 
if  the  original  problem  was  a  generalized  transportation  problem  with 


uncapacitated  arcs,  all  admissible,  then  the  check  on  total  supply  and 
total  demand  of  the  transformed  problem  provides  a  sufficient  as  well  as 
necessary  criterion  for  determining  whether  the  original  problem  has  a 
feasible  solution.  )  Finally,  our  results  imply  that  previous  theoretical 
developments  for  pure  networks  can  be  applied  to  generalized  networks  of 
suitable  rank. 

1.  I  Problem  Formulation 

Rigorous  definitions  of  networks  have  appeared  many  times  in  the 
literature,  and  we  will  not  reiterate  them  here.  Speaking  less  formally, 
a  pure  network  may  be  regarded  as  consisting  of  m  nodes  or  junction  points, 
which  are  connected  pairwise  by  a  collection  of  n  directed  arcs.  Although 
it  is  not  necessary  for  all  pairs  of  nodes  to  be  joined,  we  shall  assume  for 
convenience  that  the  network  is  connected;  i.e.  ,  there  is  a  path  consisting 
of  a  sequence  of  arcs  of  the  network  (possibly  with  some  arcs  having  their 
orientation  reversed)  which  connects  any  two  given  nodes.  Any  such  network 
may  be  concisely  described  by  its  node -arc  incidence  matrix  E.  E  is  an 

m  x  n  matrix  with  entries  E  given  as  follows: 

PH 

(-1  if  arc  q  begins  at  node  p  (node  p  is  its  tail) 

Epq  A  0  if  arc  q  does  not  touch  node  p 


►  1  if  arc  q  terminates  at  node  p  (node  p  is  il  s  head) 
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Each  column  of  E  has  two  nonzero  entries,  -1  and  +1,  iden¬ 
tifying  a  directed  arc  (i,j)  from  node  i  to  node  j.  For  each  such  admis¬ 
sible  (existing)  arc  (i,  j )  in  the  network,  we  also  define  the  following  items: 

a)  xy  is  the  flow  from  node  i  to  node  j  • 

b)  c^|  is  the  cost  of  sending  a  single  unit  of  flow  from  node  i  to 
node  j  . 

c)  b-  is  the  amount  of  supply  or  demand  at  node  i,  where  supplies 
are  denoted  by  negative  quantities  and  demands  are  denoted  by 
positive  quantities. 

Letting  N  be  the  set  of  all  admissible  arcs  (i,  j )  in  a  network,  then  a 
capacitated  pure  network  problem  can  be  stated  in  the  form: 

Problem  1 

Minimize  £  c. .  Xj. 

(i,  j)®N 

subject  to: 

-2  xii  +  E  x-.  =  b  ,  i  =  1, 2 . m 

(i,j)eN  J  (j.i)cN  J 

L.j  £  x^  £  ,  (i,  j)  e  N 

where  L- ^  and  U-  are  respectively  the  lower  and  upper  bounds 

on  the  amount  of  flow  on  the  arc  (i,  j)  [U^  need  not  be  finite]  and 

m 

where  we  assume  2j  b-  =  0. 

i=l  1 

In  contrast  to  a  pure  network,  a  generalized  network  contains  nonzero 
entries  in  its  incidence  matrix  which  are  not  necessarily  -1,  although  each 


-r>- 

column  of  the  matrix  still  contains  only  two  entries  of  opposite  sign.  We 
shall  denote  the  incidence  matrix  of  a  generalized  network  by  G.  With 
appropriate  scaling  of  each  column  of  G  and  its  associated  cost  c .  ^ ,  the 
negative  element  of  each  column  can  be  set  equal  to  -1.  Letting 
denote  the  positive  component  of  such  a  scaled  column,  the  flow  along  the 
arc  (i.j)  may  be  viewed  as  being  subjected  to  amplification  (or  attenuation) 
by  the  factor  kj  ^  (i.  e.  ,  the  amount  flowing  into  the  terminal  node  j  of  arc 
(i.j)  is  k-j  times  the  amount  that  leaves  the  initial  node  i).  The  special 
purpose  solution  procedures  [1,  2,  3,  4,  5,  6,  10,  11,  12,  13,  14]  assume 
that  this  scaling  of  the  columns  has  been  performed  before  employing  their 
solution  techniques.  The  procedure  to  be  developed  here  automatically 
performs  this  scaling  as  a  byproduct  of  its  other  operations.  For  the  case 
of  full  row  rank,  the  procedure  also  creates  a  list  of  the  arcs  which  were 
successfully  transformed  into  "pure  arcs";  i.e.,  those  whose  columns 
contain  nonzero  entries  of  1  and  -1  in  the  scaled  incidence  matrix.  Such 
a  list  can  be  used  to  computational  advantage  in  implementing  special 
purpose  solution  procedures. 

Denoting  the  nonzero  entries  of  a  column  of  G,  corresponding  to 
the  directed  arc  (i,j)  from  node  i  to  node  j,  by  -f-j  and  h^  (where  f. ^  >  0 
and  hjj  >  0),  the  capacitated  generalized  network  problem  can  be  stated  in 


the  form: 
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Problem  2 

Minimize  2  c..  Xj' 

(i,j)eN  J 

subject  to: 

-E  f .  x .  +  Tj 

(i,j)eN  13  13  (j,i)eN 


h. .  x^ 
J»  J1 


bi 


Lil  *  Xjj  s  .  <i.  j  >  «  N 

m 

where  S  b.  is  not  necessarily  equal  to  zero. 
i  =  1  1 


i  =  1, 2,  . .  .  ,  m 


Note  that  if  the  constraints  of  Problem  2  are  given  as  inequalities  rather 
than  equalities,  it  is  possible  to  transform  the  problem  into  the  format  of 
Problem  2  by  adding  absorbing  loops  (see  [11]),  although  this  may  alter 
the  rank  of  the  incidence  matrix.  In  this  connection  it  is  important  to  point 
out  that  our  procedure  can  be  applied  to  the  inequality  form  of  Problem  2 
as  well  as  to  the  equality  form.  Thus,  if  the  incidence  matrix  for  the  in¬ 
equality  form  has  rank  m-1,  one  can  still  obtain  an  equivalent  pure  network 
problem  which  can  be  converted  into  the  format  of  Problem  1  using  the 
standard  transformation. 


1. 2  Scaling  Procedure  and  Equivalence 

The  procedure  developed  in  this  section  simultaneously  determines 
the  rank  of  an  mxn  incidence  matrix  of  a  generalized  capacitated  network 
and  produces  an  equivalent  network  problem  whose  negative  entry  in  each 
column  is  -1.  If  the  original  incidence  matrix  has  rank  m-1,  this  procedure 
also  yields  a  pure  network  problem. 
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Our  procedure  can  be  viewed  as  a  labeling  method  in  which  each  row 
and  column  of  the  incidence  matrix  is  labeled  exactly  once.  The  generation 
of  these  labels  proceeds  by  scaling  one  constraint  (row  of  the  incidence 
matrix)  and  several  variables  (columns  of  the  incidence  matrix)  at  each 
iteration. 

The  scaling  procedure  creates  two  vectors  R  and  K  where  the  i^1 
component,  r.,  of  R  represents  the  scaling  factor  for  row  i  and  the  j*h 
component,  kj,  of  K  represents  the  scaling  factor  for  column  j.  Initially 
these  vectors  are  set  to  zero,  with  rj  =  0  (k  =  0)  indicating  that  row  i 
(column  j  )  has  not  been  scaled.  When  the  procedure  terminates  all  com¬ 
ponents  of  R  and  K  are  positive.  During  the  intermediate  stages  of  the 
procedure,  R  is  a  nonnegative  vector  while  the  components  of  K  may  be 
positive,  zero,  or  negative.  This  sign  convention  on  the  components  of 
K  has  the  following  interpretation: 

(a)  k|  0  indicates  that  neither  of  the  rows  containing  the  nonzero 
entries  in  column  j  have  been  scaled  (hence,  the  r^  values 
for  these  rows  are  zero). 

(b)  kj  <  0  indicates  that  exactly  one  of  the  rows  containing  the 
nonzero  entries  in  column  .J  has  been  scaled. 

(c)  kj  >  0  indicates  that  both  of  the  rows  containing  the  nonzero 
entries  in  column  j  have  been  scaled. 
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The  sealing  procedure  rnay  then  be  stated  as  follows 

1.  Set  R  =  0  anti  K  =  0.  Set  NC  0.  (At  the  termination  of  the 
scaling  procedure,  this  set  will  contain  the  indexes  of  those 
columns  whose  entries  arc  not  ±1.  ) 

2.  Scale  row  1  of  Cl  by  one  (i.e.,  set  Cl*  1  •  Cl*  and  set  rj  1). 

2.  Scale  each  column  having  a  nonzero  emry  in  row  1  by  the  reci¬ 
procal  of  its  absolute  value  in  row  1.  (That  is,  set  G,  r  — - - C 

J  l«.jl 

and  k,  =  -  — r— — r  .  for  all  j  such  that  g, .  i  0.) 

1  Ul,l  IJ 

4,  If  K  >0,  go  to  6  (by  connectedness,  all  rows  anti  columns  have 

a  , 

been  scaled).  If  K  ^  0  select  j  such  that  •  0  and  identify  > 
such  that  g.  •  *  t  0  and  r...  =  0.  Scale  row  i  by  setting 

r  .  =  __ - -  G1  =  t - i—  G*  .  and  update  K  by  setting  k, 

1  Uj  |Ri*j«l  J 


5.  Scale  each  previously  unsealed  column  with  a  nonzero  entry  in 
roiv  i  of  O  by  scttinR  k  .  1  ,  amlG.  ' - O, 

>  TcI^T  ■'  Ui.ji  J 

for  each  j  such  that  g . .  f  0  and  k  0.  Then  update  K  and 

J  J 

NC’  by  setting  kj  -kj  for  •  uch  j  such  that  g.  ^  /  0  and  k  ‘  0, 
and  setting  NC  “  NCJ{.|}  lor  each  of  these  latter  j  for  which 
g.  ^  /  1.  Then  return  to  Step  4. 


6.  If  NC"  0,  go  to  stop  7.  If  not,  for  oaoh  .)  €  NC  identify  tho 


indox  i  such  that  gj  j  <  0  and  if  gj*.  i  -1  sot 

k'  75^1“'  °>  TS^T'v 


7.  Stop. 

Onoo  tho  proooduro  torminatos  with  H  >  0  anil  K  >  0,  tho  scaling  is 
completed  by  multiplying  each  b.  by  its  corresponding  component  of  It. 
and  by  multiplying  each  cjj  and  dividing  oaoh  I.j  and  l'.  hv  its  associated 
oom|)onont  of  K.  Otherwise,  tho  problem  is  unaltered;  i.c.  ,  uneapae itated 
problems  remain  uncapacitated,  transportation  networks  remain  transpor¬ 
tation  networks,  cost  entries  retain  their  original  sign,  sources  and  sinks 
maintain  their  identity.  (This  procedure  can  be  applied  directly  to  network 
representation  rather  than  to  the  incidence  matrix.  ) 

Theorem  NC  0  when  the  sealing  procedure  terminates  if  and 
only  if  t h« ■  original  incidence  matrix  <1  has  rank  m-l. 

Proof  I'irst,  suppose  NC'  0  when  the  scaling  procedure  terminates. 
Then  every  column  of  the  scaled  incidence  matrix  contains  two  nonzero 
entries  whose  values  are  *1  and  -1  as  a  result  of  the  following  ohs<  r- 
vations  each  column  of  the  original  incidence  matrix  contained  exactly 
two  entries  ol  opixtsitc  sign,  the  scaling  weights  arc  all  positive,  and  th« 


si  ding  procedure  does  not  terminate  until  there  exists  no  column  which  has 


imi'  of  its  nonxcro  elements  in  a  sealed  row  and  the  other  in  an  unsealed 
row  (see  step  4).  !iy  connectedness,  this  implies  all  row’s  have  been  scaled 
Furthermore,  each  time  a  row  is  scaled,  each  nonzero  entry  in  the  row  that 
is  also  in  a  scaled  column  is  checked  to  see  if  it  is  II.  II  not,  NC  is  aug¬ 
mented  to  include  this  entry's  column  index.  Also,  any  nonzero  entry  in 
the  row  that  is  not  in  a  scaled  column  is  scaled  to  one  (see  step  5).  This 
demonstrates  that  NC  0  implies  that  all  nonzero  entries  of  the  scaled 
incidence  matrix  are  +1  and  -1,  which,  of  course,  means  that  the 
scaling  procedure  has  transformed  the  generalized  network  problem  into 
an  equivalent  pure  network  problem,  and  the  rank  of  Cl  is  m-1. 

Next  assume  that  the  rank  of  G  is  m-1.  Charnes  and  It;  ike  |4|  show 
that  if  any  row  of  Cl  is  deleted,  the  remaining  rows  are  linearly  independent 
Thus  if  G  has  rank  m-1,  any  row  of  G  can  be  formed  as  a  unique  linear 
combination  of  the  other  rows  in  G.  The  scaling  procedure  determines  this 
unique  combination  for  row’  1.  In  particular,  upon  applying  Steps  2  and  2  of 
the  procedure  to  scale  each  nonzero  entry  in  row  1  to  be  1 1 .  Steps  4  and  5 
are  designed  to  determine  the  unique  linear  combination  of  the  remaining 
rows  which,  when  added  to  row  1,  is  equal  to  zero.  To  accomplish  this, 
each  of  the  il  entries  currently  in  tow  1  after  applying  Step  2  must  be 
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"offset"  by  scaling  the  appropriate  row  so  that  the  unique  nonzero  entry 
which  is  contained  in  the  same  columns  as  the  -1  entry  in  row  1  is  itself 
scaled  to  -1.  This  uniquely  determines  some  r^,  and  successive  appli¬ 
cation  of  the  above  approach  yields  Steps  4  and  5  and  all  rj  values.  If 
NC  /  0  at  the  termination  of  the  scaling  procedure,  then  at  some  applica¬ 
tion  of  Step  5  there  exists  a  nonzero  entry  g^5ls-  in  row  i*  different  from 
+1  whose  associated  nonzero  column  entry  has  already  been  scaled  to  ^1. 
Moreover,  the  weight  r p;:  has  been  uniquely  identified  at  Step  4  to  create 
a  linear  combination  of  the  remaining  rows  equal  to  the  negative  of  row  1, 
if  this  is  possible,  and  the  fact  that  gi=:;j  cannot  "offset"  its  associated 
nonzero  column  entry  establishes  that  such  a  linear  combination  does  not 
exist.  Consequently,  NC  must  be  equal  to  the  null  set  if  G  has  rank  m-1. 

The  scaling  procedure  developed  in  this  section  is  quite  simple  and 
easy  to  implement.  Clearly,  as  a  result  of  the  theorem,  each  column  j 
of  the  transformed  G  contains  a  -1  and  some  positive  number  p^,  where 
p^  1  for  at  most  the  cardinality  of  the  set  NC.  An  interesting  question 
arises  as  to  how  the  cardinality  of  NC  can  be  minimized.  A  procedure  for 
accomplishing  this  would  allow  one  to  take  increased  computational  advan¬ 
tage  of  the  pure  network  structure  in  a  generalized  network  problem  that 
could  not  be  completely  reduced  to  a  pure  network.  It  is  reasonable  to 
suppose  that  "largely  pure"  networks  can  be  solved  more  efficiently  than 
those  which  are  saddled  with  a  greater  percentage  of  "generalized"  arcs. 


We  defer  such  considerations  to  another  paper. 


Consider  the  uncapacitated  generalized  network  problem: 
Minimize  3/2x^  +  ^x13  +  ®x23  +  ^x24  +  *x25  4  +  *x34 

+  2/3  +  16x^  +  5x^^ 


Subject  to 
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2  5 

x36 

4 

- 

»  • 

U 

X34 
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x45 
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x  ?0 

1.) 


The  standard  scaling  of  this  problem,  to  introduce  a  -1  entry  in  each 
column  (which  is  required,  for  example,  in  the  procedures  of  [  1,  2,  6,  1 1,  12,  13] 
is  accomplished  by  multiplying  the  columns  by  2,  1/4,  1/3,  1/4,  2,  1/2,  2, 

3,  1/4,  and  1/5,  respectively,  to  yield  the  equivalent  generalized  network 


problem : 


-  ]  3  - 


Minimize  3x12  •  lx13  +  3x23  *  *x24  +  «x36  4  2x34  +  2x4  ^  4  4x45  *  lx46 
Subject  to 


This  problem  is  graphically  depicted  in  Figure  1.  The  number  appearing 
in  the  rectangle  on  each  arc  is  the  cost  per  unit  of  flow  along  the  arc,  and 
the  number  in  the  semicircle  is  the  "multiplier  value"  for  the  are. 
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Alternative  ly,  applying  the  scaling  procedure  of  the  preceding  section 
to  the  original  incidence  matrix,  we  obtain  the  following  pure  network  problem 

Minimize  3y12  +  ly13  +  3/2  y23  +  l/2y24  +  ly25  +  4y36  +  ^34 

+  ly43  +  2y45  +  1/2y46 


Subject  to 


y  ^  0  where  R  =(1,  2,  2,  2,  8/3,  2)  and  the  components  of  K  are 
ij 

given  by  the  coefficients  of  the  y  in  the  identities  Xj9  =  2y19, 

x  1 3  =  1  ^ 4  y  1 3 '  x23=1^6y23’  x24  =  1^8  v24'  x25 = ly25* 
x36  =  1/4  y  36’  x34  =  ly34-  x43  =  ;^2y43*  x45  =  l/8y45*  and 
X46  ~  1  ;  10  y46’ 

The  fact  that  total  supply  equals  total  demand  in  this  new  representation  es¬ 
tablishes  the  possible  existence  of  a  feasible  solution  lor  this  problem. 
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To  appreciate  the  difference  in  computational  effort,  required  to 
solve  these  equivalent  statements  of  the  original  problem  the*  reader  is 
encouraged  to  apply  an  appropriate*  method  (e.  g.  ,  that  of  |  111)  to  the 
generalized  network  formulation  and  to  apply  a  standard  "minimum  cost 
flow"  algorithm  (such  as  the  out-of-kilter  method  |14|>  to  the  pure  network 
formulation. 
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